ON THE SECOND-ORDER CORRELATION FUNCTION 
OF THE CHARACTERISTIC POLYNOMIAL 
OF A HERMITIAN WIGNER MATRIX 

F. GOTZE 1 AND H. ROSTERS 

Abstract. We consider the asymptotics of the second-order correlation function 
of the characteristic polynomial of a random matrix. We show that the known 
result for a random matrix from the Gaussian Unitary Ensemble essentially 
continues to hold for a general Hermitian Wigner matrix. Our proofs rely on 
an explicit formula for the exponential generating function of the second-order 
correlation function of the characteristic polynomial. 



1. Introduction 

The characteristic polynomials of random matrices have attracted considerable 
interest in the last years, a major reason being the striking similarities between 
the (asymptotic) moments of the characteristic polynomial of a random matrix 
from the Circular Unitary Ensemble (CUE) and the (asymptotic) moments of the 
value distribution of the Riemann zeta function along its critical line (see Keating 
and Snaith [KS]). These findings have inspired several authors to investigate the 
moments and correlation functions of the characteristic polynomial also for other 
random matrix ensembles (see e.g. Brezin and Hikami [BHT] |BH2] . Mehta and 
Normand [MNJ, Strahov and Fyodorov [SF], Baik, Deift and Strahov 
[BPS] . Borodin and Strahov [BSj). 

In this paper, we consider the second-order moment and correlation function of 
the characteristic polynomial of a general (Hermitian) Wigner matrix: Let Q be 
a probability distribution on the real line such that 

JxQ(dx) = 0, a := J x 2 Q(dx) = 1/2, b := J x 4 Q(dx) < oo , (1.1) 

and let (Xu/^) i&h (X^ c ) i<jt jj eN and (Xj™ )i<j,i,jeN be independent families 
of independent random variables with distribution Q on some probability space 
(Sl,F,F). Also, let X tj := X^ c +iX}f and X j{ := Xf- e -\X}f for i < j, i,j G N. 
Then, for any N £ N, the (Hermitian) Wigner matrix of size N x iV is given by 
Xn = {Xij)i<i j<N 'j an d the second-order correlation function of the characteristic 
polynomial is given by 

f(N; ix, v) := E(det(A0v - »I N ) ■ det(X N - vI N )) , (1.2) 
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where //, v are real numbers and Ijv denotes the identity matrix of size N x N. 
We are interested in the asymptotics of the values f(N; fijy, z^v) as N — > oo, where 
/i^v? ^tv depend on A" in some suitable fashion. 

In the special case where Q is the Gaussian distribution with mean and variance 
1/2, the distribution of the random matrix Xpj is the so-called Gaussian Unitary 
Ensemble (GUE). (See e.g. Forrester [E] or Mehta |Mej and the references 
cited therein. However, let it be noted that these authors use the variance 1/4 
instead of 1/2, so that we have to do some rescalings when using their results.) 
A remarkable feature of the GUE is that the joint distribution of the eigenvalues 
of the random matrix Xn is known explicitly: It is given by 

N 

P N (dX l ,...,d\ N ) = Z N 1 ■ n (\ j -A i f-l[e- x i/ 2 A N (d\ 1 ,...,d\ N ), 

l<i<j<N 1=1 

where A N denotes the A^-dimensional Lebesgue measure on M. N and denotes 
the normalizing factor Zjq = (27t) N / 2 • Hfc=i ^ • Thus, the correlation function of 
the characteristic polynomial can be written as 

„ N N 
f GVE (N;fi,u):= / TT(Ai - fi) TT(Ai - v) P N (d\\, dX N ) , 

from which it follows (see e.g. the proof of Proposition 4.3 in Forrester [Fo]) 
that 

Here, the scalar product (•, •) is given by (<p,ip) ■= f^™(p(x)ip(x)e~ x2 ^ 2 dx, 
the pk are the monic orthogonal polynomials associated with this scalar product 
(i.e., up to scaling, the Hermite polynomials), and the kernel Kjy is given by 

K N (x, y) : = e ~(*W)/4 y P k; i{*)P k -M 

{Pk-i,Pk-l) 

Using this representation, it is possible to obtain asymptotic approximations of 
the values /gue^; MiV, v n) from the corresponding asymptotics of the Hermite 
polynomials (see e.g. Section 8.22 in Szego [Hi]). It turns out that 



tv^ooV 2N Nl J \ ' y/N* y/N J ir(n-u) 
(see e.g. the proof of Proposition 4.14 in Forrester [Fb]) and, for £ G (—2, +2), 



(see e.g. the derivation of the semi-circle law in Chapter 4.3 in Forrester [Fo] ) . 
More generally, it is known that, for £ £ (— 2, +2), 



Jim W -J L . e -^ 2 /2 . f ( N ^ + ^ + _ 

Ar-.ac V 2ttN Nl 7GUE V \/ATe(£) ^(07 

e(0 



COH-^/a^e) sin7r(/i-i/) 



vr(^ - i/) 
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(see e.g. Section 2.1 in Strahov and Fyodorov [SF]), where g(£) := ^7 \/4 - £ 2 
denotes the density of the semi-circle law. Note that this formula includes the 
preceding two formulas as special cases. 

Even more, it turns out that a similar result holds for the correlation function (of 
any even order 2, 4, 6, . . .) of the characteristic polynomial of a random matrix from 
the larger class of unitary-invariant ensembles (see e.g. Section 2.1 in Strahov and 
Fyodorov [SF]). In this respect, it is interesting to note that the emergence of 
the sine kernel is "universal" in that it is independent of the particular choice 
of the potential function of the unitary-invariant ensemble. In contrast to that, 
most of the other factors in the above result for the GUE have to be replaced by 
potential-specific factors. 

It is well-known that the GUE is a special case not only of a unitary-invariant 
ensemble but also of a (Hermitian) Wigner ensemble as described at the beginning 
of this section. The purpose of this paper is to show that the above result for 
the GUE can also be generalized in this direction. More precisely, our main result 
is as follows: 

Theorem 1.1. Let Q be a probability distribution on the real line satisfying 
let f be defined as in let £ S (—2, +2), and let |j,i/6l Then we have 



N-+oo\2ttN N\ V VNg(£) VNq(£)J 

= ex V (b-i).e^»^O. e (0 SU> 771,1 



it(h - v) 



where := ^-y4 — £ 2 and sinO/0 := 1. 

Specifically for the Gaussian distribution with mean and variance 1/2, we have 
b = |, so that we re-obtain the above result for the GUE. 

Furthermore, we see that for general Wigner matrices, the appropriately rescaled 
correlation function of the characteristic polynomial asymptotically factorizes into 
the universal sine kernel, a universal factor involving the density of the semi-circle 
law, and a non-universal factor depending only on the fourth moment b, or the 
fourth cumulant b — §, of the underlying distribution Q. 

In particular, it follows immediately that if we normalize the correlation function 
of the characteristic polynomial by means of its second moment, we obtain the 
following universality result: 

Corollary 1.2. Under the assumptions of Theorem \l.l[ we have 

j. E(Ay(g,/x) D N (£,v)) _ smTT(fi-i^) 

N-^ ^/ED N (Z,H) 2 ^EDnJ^? ~ 7r(/i-i/) ' 

where D N (£, r?) := det (Xjy - (y/NZ + 



'v 
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Moreover, it can be shown that the asymptotics remain unchanged if we replace 
the correlation function f(N; [i, v) by the "true" correlation (in the sense of prob- 
ability) of the characteristic polynomial, 

f(N; ii, v) := E((det(X JV - fiI N ) -Edet(X N - nI N )) 

■ (detpsTjv - vI N ) - EdetpQv - vIn))) ■ (1-3) 
We then have the following result: 

Proposition 1.3. LetQ be a probability distribution on the real line satisfying 
let f be defined as in tl.3\) . let £ 6 (—2, +2), and let fi, v £ R. Then we have 



lim w _L ._L . e -^ 2 /2 .j( N .yN£ + v_yNt= + ^^_) 



3\ . t A0*+»)/2e® sinn(fi-v) 



exp 6 - | • fW+wm . gg) 



TT(fl - v) 



where g(£) := ^:a/4 — ^ 2 and sin 0/0 := 1. 

Similarly as before, normalizing the correlation of the characteristic polynomial 
by means of its variance leads to the following universality result for the correlation 
coefficient: 

Corollary 1.4. Under the assumptions of Proposition HOI we have 
nm E(D N (£,iJ,)D N (Z,v)) _ smir(ti-v) 

where D N (£, V ) := det (x N - (VN£ + /jv)-Edet (x^-(ViVe + I N 

The proofs of the above-mentioned results on the correlation function of the 
characteristic polynomial of a random matrix from the GUE (or another unitary- 
invariant ensemble) heavily depend on the special structure of the joint distribution 
of the eigenvalues. However, such a structure seems not to be available for general 
Wigner matrices. 

Instead, we start from recursive equations for the correlation function of the 
characteristic polynomial (as well as some closely related correlation functions), 
derive an explicit expression for the associated exponential generating function 
and deduce all our asymptotic results from this expression. For the determinant of 
a real symmetric Wigner matrix, a similar analysis was carried out by Zhurbenko 
» 

The crucial step in our analysis is to obtain an expression in closed form for 
the exponential generating function of the correlation function. Unfortunately, 
this approach has proven successful so far only for the second-order correlation 
function of the characteristic polynomial, which explains why we do not have any 
results for the higher-order correlation functions of the characteristic polynomial. 
Note however that for any distribution Q with the first 2k moments identical to 
the Gaussian moments, the correlation function of order k of the characteristic 
polynomial must be the same as that for the GUE. 
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This paper is organized as follows. In Section 2, we start with the analysis of 
the recursive equations for the correlation function of the characteristic polynomial 
and derive the explicit expression for its exponential generating function. Sections 
3 and 4 are devoted to the proofs of Theorem 1 1 . 1 1 and Proposition [L3l respectively. 

Throughout this paper, K denotes an absolute constant which may change from 
one occurrence to the next. 

Acknowledgement. We thank Mikhail Gordin for bringing the connection 
between the correlation function of the characteristic polynomial of the GUE and 
the sine kernel to our attention. 



To simplify the notation, we adopt the following conventions: The determinant 
of the "empty" (i.e., 0x0) matrix is taken to be 1. If A is an n x n matrix and 
z is a real or complex number, we set A — z := A — zl n , where I n denotes the 
n x n identity matrix. Furthermore, if A is an n x n matrix and i\, . . . ,i m 

and 

ill • • • ; jm are families of pairwise different indices from the set {1, . . . , n}, we write 
j ^[ii,...,i m :ji,...,j m ] f or tfig [n — m) x (n — m)-matrix obtained from A by removing 
the rows indexed by i\, . . . , i m and the columns indexed by j%, . . . ,j m . Thus, for 
any n x n matrix A = (d«)i<i ,< n (n > 1), we have the identity 



det(^) = a i>n a n>j det(A^ n ^) + a n , n det(^ :n l) , (2.1) 



as follows by expanding the determinant about the last row and the last column. 
(For n = 1, note that the big sum vanishes.) 

Recall that we write X/v for the random matrix (Xjj)i<jj<7v, where the are 
the random variables introduced below (jl.ip . We will analyze the function 



2. Generating Functions 



n-l 



f(N; pi, v) := E (detpfTv - fj.) ■ det{X N - v)) 



{N > 0) . 



To this purpose, we will also need the auxiliary functions 



A x i(iV;^) 

foi(N;n,v) 



E(det((X N - / u)I 1:1 l) • det((X N - u)^)) 
E(det((X7v - fj,)^) ■ det((X N - v f^)) 
E(det(XAT_i - fj,) ■ det(X N - v)) 
E(det(XAr - n) ■ det(Xiv_i - v)) 



(N>2), 

(N>2), 

(N>1), 
(N > 1) . 



Since [i and v can be regarded as constants for the purposes of this section, we will 
only write f(N) instead of f(N; (i, v) in the sequel, etc. 
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We have the following recursive equations: 
Lemma 2.1. 



/(0) = 1, 

f(N) = (1 + (mv) f(N - 1) + (26 + |)(iV - 1) f(N - 2) 
+ (iV-l)(JV-2)/u(JV-l) 
+ (N-l)(N-2)f*{N-l) 
+ v(N - 1) f 10 (N - 1) 

+ fi(N-l)f 01 (N-l) (N>1), (2.2) 

fn(N) = fivf(N - 2) + (N - 2) f(N - 3) 
+ (AT_2)(Ar_3)/ 11 (Ar_2) 
+ v{N -2)f w {N -2) 

+ /i(JV-2)/oi(JV-2) (JV>2), (2.3) 

/*(7V) = /(iV - 2) + (iV - 2) /(iV - 3) 

+ (JV — 2)(iV — 3) fii(N — 2) (N>2), (2.4) 

f 10 (N) = -(N-l)f 01 (N-l)-uf(N-l) (N>1), (2-5) 

f m (N) = -(N-l)f 10 (N-l)- f if(N-l) (N>1). (2.6) 

For the sake of clarity, note that these recursive equations may contain some 
terms which have not been defined (such as fn(N — 1) for N = 1), but this is not 
a problem since these terms occur in combination with the factor zero only. 

Proof. We begin with the proof of (|2.2p . For N = 0, the result is clear. For N > 1, 
we expand the determinants of the matrices (Xjy — fi) and (Xn — v) as in (12. lj) 
and use the independence of the random variables Xy = Xji (i < j), to the effect 
that 

/(AO 
N-1 N-1 

= X X {-^) l+j+k+l E (X iiN X NJ X ktN X Ntl ) ■ E (det(X Ar _ 1 - • det(X Ar _ 1 - 
i,j=l fc,Z=l 
AT-1 

+ 2 E (X itN X Nd ) • E (Xiv,iv - i/) • E (detpOv-i - • det(Xv-i 

IV-1 

+ XI ("l) fe+Z+1 E (X KN X N)l ) ■ E (Xiv,iv - fi) ■ E (det(Xv-i - n) ■ det(Xv-i - i 

k,l=l 

+ E ((X NtN - n)(X NsN — v)) • E (det(Xjv-i - jx) ■ det(Xv-i - u)) . 



Since the (complex-valued) random variables Xjj = X/, (i < j) are independent 
Uj) = (» < j) and E(X§) 



with E(Xjj) = (f < j) and E(X 4 2 -) = (i < j), several of the expectations vanish, 



Characteristic Polynomials of Hermitian Wigner Matrices 7 

and the sum reduces to 
f(N) 

= (EX^jy + [iv) ■ E (det(Xjv_i - /i) • det(X iV -i - i/)) 

+ ^ E|X 4 , iV | 4 -E(det(X i v_ 1 -/i)^-det(X i v_i-^) [fc:/] ) 

+ ]T E|X^| 2 -E|X fciJ v| 2 -E(det(X i v_ 1 - M )[^-det(X^ 1 -^)^]) 
i=j^k=l 

+ ]T E|X i , i v| 2 -E|X fc ^| 2 -E(det(X i v_ 1 - M )[^-det(X w _ 1 -^)[ fe: ']) 

i=l^=j=k 

+ v^2E\X iiN \ 2 • E (detpOv-i - ^) [ ^' ] • det(Xiv-i - vj) 
i=j 

+ /^E|X fcj7V | 2 ■ E (det(Xjv_i - m) • det(Xiv_i - ^) [fe: ' ] ) . 

fc=Z 

(P now follows by noting that EX%^ N = 1, E\X iyN \ 2 = 1, E\X i)N \ 4 = 26 + 1, and 
by using symmetry. 

To prove (j2.3l) . we apply the analogue of (|2.1|) for the first row and the first 
column to the matrices (Xjv — [ip 1 * and (X* — ^)' 2:2 ^. Using similar arguments 
as in the proof of (|2.2D afterwards, we obtain 

/n(iV) 

N N 

= E E E (X i)2 X 2j ) • E (X M X y ) • E (det(Xjv - ,,) '- 2 - 1 - 2 -.' • dettX* - y)lW> 2 >^ 

i,j=3 k,l=3 
N 

+ (-iy +j+1 E (X i)2 X 2J ) • E (X 1A -v)-E (det(X N - . det (X N - u)^ a ^j 

i,j=3 
N 

+ E ("l) fc+m E (X fe)1 X M ) • E (X 2 , 2 - /i) • E (de^X* - ^) [1 ' 2:1 ' a] ■ det(X* - v^Mtf) 

k,l=3 

+ E (X 2>2 -fi)-E (X 1;1 — u)-E (det(X N - n) [1 ' 2ll ' 2] ■ det(X iV - j,)! 1 ' 211 ' 2 ]) 
= ^ • E (det^ - m) [1 ' 2:1,2] • det(X w - v )lW>A 

+ E \ X ^\ 2 ■ E I X mI 2 • E (det(Xjv - • det(X w - vf^M^ 

i=j=k=l 

+ E \ X ^\ 2 ■ E I X mI 2 • E (det(XAr - n) [1 ' 2M ' 2 ' j] ■ det(X N - yf^M^ 

i=j^k=l 

+ ^E|X i>2 | 2 • E (det(X N - Ai )[i.2.fci,2j] . det( x N - ^[1,2:1,2] j 
+ ^J>|X M | 2 • E (detOX* - m) [1,2:1,2] • det(X* - 

k=l 

and hence ()2.3|) . 
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The proof of (|2.4p is similar to that of (|2.3p : 

N N 

- E E (-lT +J+k+l E (Xi.iXi,,) • E (X 2ji X fc , 2 ) • E (det(Xjv - ^I 1 ' 2 * 1 ^ • det(X* - u f>W<^ 

i,j=3 k,l=3 
N 



+ E E (X^xXsj) • E (X 1)2 ) • E (det(Xiv - ^ ia ^ • det(X,v - v) [1 ^ 2] ) 

+ E ("l) fc+m E (Xfc.aXi,,) • E {X 2 ,i) ■ E (det(Xjv - /i) [1 ' 2:1 ' 2] • det(Xjv - ^A*:!,^ 



k,l=3 

+ E (X 2 ,iXi j2 ) • E (det(Xjv - m) [1 ' 2:1,2] • det(X7v - ^) [1,2:1 ' 2] ) 
= E (det(^v - m) [1,2:1,2] • det(Xjv - uf^) 

+ £ EIX^I 2 • E|X fc , 2 | 2 • E (det^-i - pfWM • det(X^! - yf^M 

i=l=k=j 

+ J] E|X i5l | 2 • E|X fc , 2 | 2 • E (det(X N ^t - ^[W^j] . det(^_! - ^AfciAq 

For the proof of (|2.5p . we expand the determinant of the matrix (Xn — v) as 
in (|2,ip and use similar arguments as above to obtain 

N-l 

fxo{N) = E E (X k , N X Nt i) ■ E (det(X iV -i - /u) • det^-i - 

fc,/=i 

+ E (X N>N - v) ■ E (det(XAT_i - //) • det(X7v_i - 1/)) 

= - E E|X fcj7V | 2 • E (det(Xjv-i - fi) • det(Xjv_i - 
k=i 

- v ■ E (det(Xjv_i - /i) • det(Xjv-i - 1/)) 



and hence (12.51) . 

The proof of ()2.6[) is completely analogous to that of (|2.5p . □ 

The interesting (although apparently rather special) phenomenon is that the 
preceding recursions can be combined into a single recursion involving only the 
values f(N). To shorten the notation, we put 



<N) := ^ (N> 0) 



and 



s (iV) := E ^ - fc ) > 0) • 

k=0,...,N 
k even 

We then have the following result: 
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Lemma 2.2. The values c(N) satisfy the recursive equation 

c(0) = l, (2.7) 
Nc(N) = c(N-l) + N- c(N - 2) 

+ \tv ■ (s(N - 1) + s(N - 3)) 
- (p 2 + v 2 ) ■ s(N - 2) 

+ (26 - |) • (c(N - 2) - c(N - 4)) (N > 1) , (2.8) 

where all terms c( • ) and s( • ) a negative argument are taken to be zero. 

Proof. It is immediate from Lemma |2. II that 

/(AO = / n (iV + 1) + f*(N + 1) + (26 - |)(7V - 1) • f(N - 2) 

for all N > 1. Using this relation for N — 2 instead of iV (where now N > 3), 
we can substitute fn(N — 1) + fu(N — 1) on the right-hand side of (|2.2p to obtain 

/(AO = (1 + H /(JV - 1) + (26 + i)(iV - 1) f(N - 2) 

+ (at _ i)(jv - 2) (/(AT - 2) - (26 - \){N - 3)f(N - 4)) 

+ u(N - 1) f 10 (N - 1) 
+ /i(iV-l)/ 01 (iV-l) 
= (1 + iiv) f(N - 1) + AT(A^ - 1) f(N - 2) 
+ i/(JV-l)/i (iV-l) 
+ ^(JV-l)/oi(JV-l) 

+ (26 - |) • ((AT - 1) /(AT - 2) - (N - 1)(JV - 2)(AT - 3) /(AT - 4)) 
for all N >3. Dividing by (N — 1)!, it follows that 

Nc(N) = (1 + /«/) • c(N - 1) + A^c(A^ - 2) 
+ ^•/ 10 (Ar-l)/(AT-2)! 
+ iM.f 01 (N-l)/(N-2)\ 

+ (26 - |) • (c(N - 2) - c(N - 4)) 

for all A^ > 3. (For N = 3, note that the second term in the large bracket vanishes.) 
A straightforward induction using (|2.5p and (|2.6|) shows that 

f w (N - 1) / (N - 2)! = -us(N - 2) + M s(AT - 3) , 

f 01 (N - 1) / (N - 2)! = - MS (A^ _ 2 ) + us(N - 3) , 

for all N >3, which yields the assertion for A^ > 3. 

For A^ < 3, the assertion is verified by direct calculation, also making use of 
Lemma 12.11 

c(0) = /(0) = l. 
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ic(i) = /(i) = (i + H/(o) 

= (l + /ii/)c(0) 
= c(0) + (ivs(0) . 

2c(2) = /(2) = (1 + + (26 + i)/(0) + K""/(0)) + M(-M/(0)) 

= (1 + + (26 + 2)/(0) " (M 2 + ^)/(0) 

= (1 + m^)c(I) + (26 + §)c(0) - (/i 2 + is 2 )c(0) 
= c(l) + 2c(0) + ^c(l) - (/u 2 + v 2 )c{0) + (26 - |)c(0) 
= c(l) + 2c(0) + ^s(l) - (/i 2 + z^ 2 )s(0) + (26 - |)c(0) . 



□ 



We can now determine the exponential generating function of the sequence 
(/(iV))jv> : 

Lemma 2.3. The exponential generating function F(x) := ^5v=o f(N) % / N\ of 
the sequence (f(N))i\r> is given by 



exp [av ■ - \(n 2 + v 2 ) ■ + b*x 2 

F W = (l_ x )3/2.( 1 + x )l/2 ' 

where b* := b — § . 

Proof. It is straightforward to obtain F{x) starting from (|2.7p and (12. 8p and 
using the basic properties of generating functions. For the sake of completeness, 
we provide a detailed proof. 

To begin with, recall that f(N) /N\ = c(N). Multiplying (j2T8j) by x N ~ x , summing 
over N and recalling our convention concerning negative arguments, we have 

oo oo oo 

)x N ~ 1 



N=l N=l N=2 

/ oo 



^ Nc(N)x N - 1 = ^2c(N- ljz*" 1 + Nc ( N ~ 2 

N=2 

oo oc 

£ s(N - l)x N - x + E s(N-3) 

N=l 

oo 

+ ^) £ s(A r _ 2 ) 

V=2 

/ oo 

26* ij2 c ( N - 2 ) xN ~ l - Yl c ( N - 4 ) 



^v=i V=3 / 

oo 
7V=2 

oo oo x 

AT-1 



N=2 N=4 



whence 

F'(x) = F{x) + (2xF(x) + x 2 F'(:e)) 

+ liV \±^F{x) - (/i 2 + u 2 )-^F{x) + 26* (xF(i) - x 3 F(x)) . 
1 — x z 1 — x z 

We therefore obtain the differential equation 

n.) = (ct + ""(^ - <" 2 + (T^f + 2t **) F| *> ■ 
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which has the solution 

Here, Fq denotes a multiplicative constant which is determined by ()2.7|) : 

F = exp(I(/i 2 + z. 2 )). 

We therefore obtain 

F{X) = (l-x)3/Ml + x)V 2 exp (^13^ " ^ + ^13? + 6V ) ' 
which completes the proof. □ 



3. The Proof of Theorem 11.11 
To prove Theorem ll.!) we will establish the following slightly more general result: 

Proposition 3.1. LetQ be a probability distribution on the real line satisfying U.l\) , 
let f be defined as in Kl.Hty . let {£,N~)N<m be a sequence of real numbers such that 
lim^v^oo £n/s/~N = £ for some £ G (—2, +2), and let ijeC. Then we have 

I 2tT 1 f TJ 



where sin 0/0 := 1. 

It is easy to deduce Theorem 11.11 from Proposition l3.lt 

Proof of Theorem Taking 

nr?+ 7r(u + v) , 7r(u — I/) 

Cn ■■= VN£ + — ^- - L = and 77 := 7 7 

in Proposition 13-H we have 
^ 1 



J im V 3 • TFT • exp -^V2 - + ^)/\/4-e 2 

TV^oo V JV JV! V 



exp • 



sin7r(^ — v) 



7T(/Z - I/) 

Multiplying by ^ exp ( 7r£(// + v)/ yA — £ 2 ] yields Theorem I l.li □ 
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It therefore remains to prove Proposition 13, It 
Proof of Proposition \3.1l By Lemma |2, 3\ we have 

" W;n,v) „_ exp((.fri ? -5(/' 2 + i' 2 )- 1^+1^) 



E 



-z 



\ . : (| ^ ^ (1-^.(1 + ^)1/2 

Thus, by Cauchy's formula, we have the integral representation 

lz-,,2 i ,,2\ z 2 i z,*.v2 



f(N;(M,v) 1 /-«q)^-T^T-3(^ + ^)-iZF + &VJ dz 



AH 2vri L (1 - z )3/2 . (l + z )i/2 



N+l 



(3.1) 



where 7 = 7at denotes the counterclockwise circle of radius R = Rn = 1 — 1/iV 
around the origin. (We may and do assume that N > 2 for the rest of the proof.) 
Setting jjL = £n + rj/y/N and v = £jy — ^//V^V, we have 

/* 2 1/2 2\ ^ i*2 

exp //f • ^ - o(M +v*)- - ^ + b z 

\ 1 — z A 1 — z z 



- * „2 



exp ( - r/ViV) • ^ - + r^/iV) ■ — 5 + 



1 - V ' JV ' ' ' 1 - 



' 2 / Art Z 1 l*„2 



exp ( e N ■ — z ~ (VVN) - — + b*z 



1-2 ,0 1 .* 2 



exp (±& + V 2 /N ) . exp / _if N . __ _ (r? 2 /iV) . — + b*z 



We therefore obtain 



7 ( jV; ^ + 71' ^ " V^v) = 6XP ^ + V2,N ^ 

1 /-ex P (-|^-fef-(WiV)-T^+^ 2 ) «fa 



27ri7 7 (1 - je)3/2 • (1 + 2)1/2 Z JV+1 



(3.2) 



The idea is that the main contribution to the integral in (|3.2|) comes from a small 
neighborhood of z = 1, where the function 

_ ex P (-^-ig-(WiV)-T 1 l+fe^ 2 ) 

(1 - *)3/2 • (1 + Z )V2 

can be well approximated by the simpler function 

exp(6*) ex P ("3^ ■(!-«)- (WAO • T 1 ! 



M*0 := 



V2 (1 - z) 3 / 2 

We therefore rewrite the integral in f|3.2|) as 



2vrf L V ; 



Z AT+1 



h+h + h- h (3.3) 
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with 



h-7T - h ( Re ) virkw > (3-5) 



h 



2vr ii/v^v 1 " (-Re' 

M^') — 5nv- (3.7) 



2vr ii/^jv " (/>'' 

We will show that the integral Ii is the asymptotically dominant term. 
First of all, note that since £jv G M, we have 

|exp(-i^-(l-z))| =exp(-i^-Re(l-z)) <1 (3.8) 

for any z E C with Re (z) < 1. Plugging in the series expansion 

i \ ~ (-i)V i 



exp (-(W^)- T 1 I )=E 



i\N l (i-z) 1 

and using uniform convergence on the contour 7 (for fixed N > 2 and 77 € C), 
we obtain 

Ji_ exp(6*) 1 f exp (4& • (1 - s)) dg 

ViV v/2 ^ Z! ^ +1 / 2 2vri J 1 (1 - z)*+3/2 ^+1 " 1 ' > 

We will show that for each I = 0, 1, 2, 3, . . . , 

'(-!)' ^ 1 /■ exp •(!-*)) dz 



lim , , 

TV— >oo \l\ N l+1 / 2 2m J (l- Z )l+3/2 



N+l 



z 



(2Z + 1)! 1 ? J 1 j 



To begin with, 

1 /" exp (-i^ . (x _ z )) d2 



2ni L (1 - z)'+ 3 / 2 



N+l 



Z 



2^/(i-i/iv)-ioo (l-z)<+ 3 / 2 



(3-11) 



In fact, for any i?' > 1, we can replace the contour 7 by the contour <5 which consists 
of the line segment between the points R — i \J (R') 2 — R 2 and R + i y / {W) 2 ^R 2 , 
and the arc of radius i?' around the origin to the left of this line segment (see Figure [1]). 
Now, it is easy to see that the integral along this arc is bounded above by 

1,1 1 
— • 2-kR' ■ 



2tt (i?'-iy+ 3 / 2 (R') N+1 

and therefore tends to zero as R' — * 00, whence f)3.11|) . 
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Figure 1. The contour 5. 



Next, performing a change of variables, we find that the right-hand side in (|3.11j) 
is equal to 



N l + i/2 J_ r +00 ew{-l(e N /N)-(l-iu)) Ha 
' 2vr 



(1 - m)'+ 3 / 2 (1 - ^) N+1 



Since lirn/v-KJo £,n/v N = it follows by the dominated convergence theorem that 



Jim J_ /' +00ex P(-i(^/ Ar )-( 1 -™)) 



(1 - iu) l +^ 2 (1 - ^) N + 1 

1 /•+- exp ((1 - \i 2 ) ■ (1 - iu)) 



2tt 



(1 - in)'+ 3 / 2 



du , 



which is equal to 



(i - k 2 ) m/2 i n 



T(/ + 3/2) y^F (2i + l)! 



2^+1/2 



by Laplace inversion (see e.g. Chapter 24 in Doetsch [Do]) and the functional 
equation of the Gamma function. This proves (|3.10p . 
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Let e > denote a constant such that cos t < 1 - e 2 t 2 for -7T < i < +vr. Then, 
for any a > 1, we have the estimate 



+7T 



1 



II - Re u \ 



dt 



1 



< 



_ w (1 + R 2 -2Rcost) a / 2 

+7T ^ 



_ w ((1 - R) 2 + e 2 t 2 ) a / 2 



dt 



N c 



1 



_ w (l + N 2 e 2 t 2 ) a / 2 

Neir i 
a-1 ' 1 



dt 



< fC/V Q_1 ( 1 + / — 



< KN a ~ l I 1 + 



a — 1 



(3.12) 



where IT denotes some absolute constant which may change from line to line. 
We therefore obtain the bound 



E 

1=0 



(-l)V 1 /-exp(4& •(!-*)) 



Z! JV'+V2 27TZ 



(1 - z)'+3/ 2 Z 



iV+1 



< 



< 



E 



21 



1 



1 



1 



11 N l+1 l 2 2vr /_ , |l-i?e i ^+ 3 / 2 lite**!* 



*E 



Z! 



21 / i 
1 + 



2 + 1/2 



< oo . 



uniformly in TV" > 2. Thus, the term- by-term convergence established in (|3.10p 
entails the convergence of the complete series in (|3,9p . and we obtain 



lim * 



(21 + 1)! 



2/ 



— • exp (b - f ) • 



z=o 

sin(y / 4 — £ 2 • 77) 



(v^F-t?) 



Hence, in view of (j3.2[) and (|3.3p . the proof of Proposition 13.11 will be complete 
once we have shown that the integrals I2, I3, I a are asymptotically negligible in 
the sense that they are of order o(yN). 
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For the integral we use the estimates (R = 1 — 1/N, t £ 



exp 



r2 1-Re u \ 
1 + Re u J 



l-ft 2 



exp 



lf2 



exp -(WiV) 



1 



1 - i?e i4 



< exp i^\v\ 2 /N) 

< exp ( (\ri\ 2 /N) 



1 + R 2 + 2Rcost 
1 



<1, (3.13) 



|1 - Re u \ 
1 



1 - R 



exp(|r/p 



exp [b*{Re lt Y) <exp M^Hiie'T <exp(|6*|), 



|1 + Re u \ > 1 + Rcost > 1 for cosi>0, 

|1 - i?e**| > 1 - Ecost > 1 for cos£<0. 
Using these estimates, it follows that 

|n< J_ f /2 exp(H 2 + |b*|) dt 

| _1_ f^' 2 exp(M 2 + |6*|) dt_ 
+ 2ttJ w/2 \l + Re it \ 1 / 2 R N 

1_ /-a^-V^ exp(|77| 2 + |6*|) dt_ 

+ ^ L/2 ll-^p/2 fiJV • 

Similarly as in (|3.12p . we have the estimates 



/ 1/V W|1-^|3/ 2 
r 3ir/2 



dt < KN 1/2 



U u3/2 



du < KN 1/A , 



/ 1 ^7 * < KN~ 1/2 1 + / —777 dw 

A/2 l+i^V2 " ^ U V2 y 



(3.14) 

(3.15) 

(3.16) 
(3.17) 



<K 



r 27r-l/ v / F ] 
/ 3V 2 |1-^'|3/ 2 



dt < KN 1/2 



Ne «' 



1 

372 



< iTiV 1/4 . 



(Recall our convention that the constant K may change from one occurrence to 
the next.) It follows that 

\h\ <KeM\v\ 2 + W\)N l ^ = o{sfN). 

For the integral ^3, we write 



h(z) - ho(z) 



exp [-\e N • (1 - *Q - (V 2 /N) ■ ^ 



(1 - zfl 2 



•UW -Mi)) , 
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Hz) 



exp 



4?JV i +z +0 Z 



(l + z) 1 / 2 



2?^ 1+2 + 4?AT (I+2) 2 + ° Z 



1/2 



where 

so that 
h'{z) - 

Let 

Z:={zGC|z = re**, 1 - 1/JV < r < 1, \(p\ < 1/Vn} 
and note that for z £ Z, we have the estimate 



(l + z) 1 / 2 



(1 + Z)3/ 2 

• exp 



1*2 



1 + z 



Re 



l + z 



< 4/JV. 



(3.18) 



Indeed, with z = re J</ , a simple calculation yields 

/ (1 — z) 2 \ — 1 + r cos </J + 2r 2 — r 2 cos 2<£> — r 3 cos </? 

Re = ^ , 

y 1 + z J 1 + r z + 2r cos </> 

where the denominator is clearly larger than 1 and the numerator is bounded 
above by 

r cos ip — r 3 cos ip + r 2 — r 2 cos 2</? < r(l — r 2 ) cos ip + r 2 — r 2 (l — 2y? 2 ) < 4/iV . 

Since for z = Re lt with \t\ < l/y/N, the line segment between the points z and 1 
is contained in the set Z, it follows that 

'h'(0 



h{z)-h(l) <\z-l\ sup h'{(l-a)z + a) 

a6[0;l] 



< \z — 1\ sup 



<K\z-l\ sup |^|l-C| + |b*| + lJ exp ^ 

< K\z - 1| + |6*| + l) exp (&/N + |6*|) 

<K(b*,C)VN\z-l\, 



(i-C) s 
i + C 



where the last step uses that liniAr^oo £jv/v N = £, and K(&*,£;*) denotes some 
constant depending only on 6* and ^* := (£Ar)jveN- Using (|3.8p . (|3.14p as well as 
a similar estimate as in (I3.12p . we therefore obtain 

Li/^v |l-i?e^| 3 / 2 



< K(b*,C,v) VN 

<K( b *,e,v)(i + l ^ 

< K(b*,C,v)N 1/4 = o(VW), 



h(Re u ) - h(l) 
1 



dt 
W 



_ 1/v ^v |l-i?e«|V2 
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where K(b*,^*,rj) denotes some constant which depends only on b* , £*, and r\ 
(and which may change from line to line as usual). 

For the integral I4, we can use (13.8p . (13.14j) as well as a similar estimate as 
in (|3,12p to obtain 

, J_ /■ 2 "- 1 A / ^ exp(|7 ? | 2 + |ft*|) dt_ 

( pit 1 r2n-l/VN i \ 

< i.expd,, 2 + |.* |) ^ f[ _^ * + jf f[ _^ *j 

< ^expdr/l 2 + |6*|) ( iV 1 / 2 /°° -I? du + N ' /2 H ^2 dn J 

< KeM\v\ 2 + \b*\) N 1 ^ = o(V7V) . 

This completes the proof of Proposition 13.11 □ 



4. The Proof of Proposition 11.31 

Let f(N; fj,, v) and f(N; fi, v) be defined as in (|1.2jl and (|1.3p . respectively, and 
note that 

7(^;/x,z/) 

= E((det(X7v - n) -Edet(X N - fi)) • (det(Xjv - v) - EdetpGv - v))) 

= E( det(Xiv - m) • det(XAT - 1/)) - Edetpfjv - //) • Edet(Xjv - v) 

= f(N;^u)-g(N;fi)g(N;u), (4.1) 

where 

g(N-fi) :=Edei(X N - 

for any /i£l. We will deduce Proposition 11.31 from Theorem 11.11 by showing that 
f(N; fj,, v) is asymptotically much larger than g(N; (j) g(N; u). 

To this end, we need some more information about the values g(N; [/,). Similarly 
as in Section 2, we have the following recursive equation: 

Lemma 4.1. 

0(O;/i) = l, g(N;tx) = -ng(N-l ] »)-(N-l)g(N-2;v) (N > 1) . 



Proof. For N = 0, the claim follows from our convention that the determinant 
of the empty matrix is 1. For N > 1, we expand the determinant of the matrix 
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(X/v — (jl) as in (|2.ip and use independence and symmetry to get 
N-i 

g(N; fi) = E (X itN X Nd ) ■ E (det^-i - ^ 

+ E (X NN - fj) ■ E (det(Xjv-i - fj)) 



E|X i)JV | 2 • E (det(X2v-i - fif- 



i=3 

+ E (X^vtv — /i) • E (det^jv-i - /i)) 
= -(JV - l) 5 (iV - 2; - M5 (iV - 1; M ) ■ 
This completes the proof of Lemma 14.11 □ 

It follows from Lemma |4. II that the polynomials g(N; //) coincide, up to scaling, 
with the Hermite polynomials H^{x) (see e.g. Section 5.5 in Szego |Szj). which 
satisfy the recursive equation 

H (x) = 1, H N (x) = 2xH N _ l {x) - 2{N - l)H N . 2 {x) {N > 1) . 

(Specifically for the GUE, this is well-known, see e.g. Chapter 4 in Forrester [Fo].) 
The precise relationship is as follows: 

Lemma 4.2. For any N = 0, 1, 2, 3, ... , 

g{N^) = {-l) N 2- N ' 2 H N ^/V2). 

Proof. This follows from the recursive equations for g(N; fj,) and Hn(x) by a straight- 
forward induction oniV. □ 

Due to Lemma 14. 2\ it is easy to obtain the asymptotics of the values 
g(N; yN£ + fi/y/N g(£)) from the corresponding asymptotics of the Hermite poly- 
nomials (see e.g. Section 8.22 in Szego |Szj). For our purposes, the following 
estimate will be sufficient: 



Lemma 4.3. For £ G (-2, +2), /j£R fixed, 



e- N ^/ 4 g N;VNH + 



where K(£,fM) is some constant depending only on £ and \x. 



Proof. By Theorem 8.22.9 (a) in Szego |Sz| . we have, for x = \/2N + 1 cosy?, 
e- x2 ' 2 H N {x) = 2W 2 )+( 1 /4) (A r,)i/2 (A r^)-i/4 (sin ^-i/2 

• (sin {({2N + l)/4) • (sin 2^-2^) + 3^/4) + 0{N- 1 )) , 

where the O-bound holds uniformly in <p E [e, it — e], for any e > 0. 

Combining this result with Lemma 14.21 we obtain, for N sufficiently large, 

= (-l) N 2 l l 4 (M) 1 / 2 (Ntt)- 1 /* (simp)- 1 / 2 

• (sin(((2iV + l)/4) • (sin2v?-2v?) +3vr/4) +C(iY" 1 )) , 
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where 

if = ifj\[ := arccos 

is contained in an interval of the form [— e, ir — e] with e > 0. From this, Lemma [4.31 
easily follows. □ 

After these preparations we can turn to the proof of Proposition 11.31 

Proof of Proposition \1.3i By Equation (|4.ip and Lemma [4.3l the difference between 
the left-hand sides in Theorem 1 1 . 1 1 and Proposition 11.31 is bounded by 



where K(£, fx, v) is some constant depending only on £, \i and v. Thus, Propo- 
sition [L3] follows from Theorem II .li □ 
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